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Abstract. 

^ , We consider a model for 2D electrons in a very strong magnetic field (i.e. projected 

onto a single Landau level) and a random potential V. The computation of the averaged 
Green function for this system reduces to calculating the averaged density of states. 
■ We have constructed a computer algebra program which automatically generates a 

, perturbation expansion in V for these quantities. This is equivalent to computing 

Qv^ I moments of the density of states. When is a sum of Gaussians from Poisson 

distributed impurities, each term in the perturbation expansion can be evaluated 
automatically. We have done so up to 12th order. The resulting information can 
be used to reconstruct the density of states to good precision. 



o 



^ ' 1. Introduction 



The Quantum Hall effects has continued to furnish theoretical physics with challenges 
and interesting problems for almost twenty years. Is is remarkable that the combination 
of basic physics (classical electrostatics and magnetism, and non-relativistic quantum 
mechanics) which has been known and researched upon for so long can lead to so many 
unexpected and exotic phenomena. 

In the work reported on in this paper we have been inspired by the integer Quantum 
Hall effect Ijl]] to initiate a somewhat prosaic (but we think nevertheless useful) project — 
the brute force high order perturbation expansion of impurity averaged quantities in the 
Quantum Hall system. The problem we consider is a model of non-interacting electrons 
confined to a two-dimensional layer in a very strong perpendicular magnetic field and 
a random impurity potential V. All information about the physical properties of this 
model is encoded in its various Green functions. In the absence of electron-electron 
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interactions it suffices to consider the one-particle Green function, e.g. the solution to 
an equation like 

-f^G{f,p;t)=HGif,p;t), 

where H is the one-particle Hamiltonian. Thus, G is a random quantity (being a 
functional of V) for which we can only expect to be able to compute various impurity 
averages like Ylk=i C^i^k, Pk] tk)- In this paper we shall only consider the case of i^' = 1. 
To investigate transport properties one must also consider the case of i^' = 2. 

Our perturbation expansion consists of expanding G into a series in powers of V 
(which for our model is equivalent to an expansion in powers of time), and performing 
the impurity average of each term in this series. The result of the averaging process can 
be represented by a sum of (Feynman) diagrams. The number of diagrams grows quite 
rapidly with the perturbative order, to 12'th order it reaches a few hundred thousand 
(the exact numbers are listed in Table 1). We have written a computer algebra program 
which automatically generates all diagrams to a given order. 

The use of the generated diagrams is not restricted to the Quantum Hall system. 
However, the great simplification which occurs for this system (under conditions 
explained below) is that each diagram reduces to a low-dimensional Gaussian integral 
(or the integral of a Gaussian multiplied by a polynomial), and thus can be evaluated 
analytically by computer algebra. The conditions are i) that the system is projected 
onto a single Landau level (most simply the lowest one), and ii) that all correlators 
V{ri) ■ ■ ■V{fk), can be written as a sum of Gaussians (or Gaussians multiplied by 
polynomials). The latter restriction is fulfilled if we assume V to arise from a density 
of Poisson distributed impurities with a Gaussian impurity potential. In the ^-function 
limit, when the range of the Gaussian goes to zero, this model belongs to a class of 
systems which was solved exactly for the density of states by Brezin et. al. P] (a simpler 
system, which corresponds to the additional limit of taking the impurity density to 
infinity, was solved earlier by Wegner|^. These solutions provide useful checks of our 
results and methodology. 

Given the beginning of a perturbation series, even to high orders, it is not obvious 
how one should extract the correct physical information from it, if its expansion 
parameter fail to be small. It is not even obvious that computing a few more terms 
in the perturbation expansion will be helpful. This is a problem which arises quite 
often, and for which many methods of massaging the perturbation expansions have 
been devised. Our model provide a non-trivial example of this problem, which can be 
investigated to rather high order of perturbation theory. We have considered the most 
common methods of summing infinite subsets of the perturbation series. The results 
were not encouraging. This is not surprising, in view of the fact that such subsets 
are usually selected more by their property of being (easily) summable than by their 



3 



physical importance. 

A more fruitful approach to our problem is to utilize the fact that the information 
contained in the n'th perturbative order is equivalent to the n'th moment (s) of a spectral 
function. For the averaged Green function this spectral function is actually the density 
of states (as function of energy). This is due to the fact that translation invariance 
forces the relation 



where V is the projection operator onto the given Landau level. 

There are, of course, infinitely many spectral functions which reproduce a given 
finite set of moments. However, for physical systems the assumptions of a certain 
degree of smoothness and simple asymptotic behaviour are usually reasonable. Then, 
the high order behaviour of the moments provide information about the tail of the 
spectral function, while their totality pins down its overall behaviour. We have applied 
this procedure with very satisfactory results. 

The rest of this paper is organized as follows: In section 2 we give a more detailed 
description of our model and the method of solution. In section 3 we describe some of the 
functions which naturally occur in the graphical expansion, the relations between them, 
and the number of graphs contributing to each perturbative order. This information 
has been used to check the correctness of our computer algebra programs. In section 4 
we present our calculated moments, for the lowest Landau level, and give some physical 
interpretation of the results. In section 5 we discuss the exact results of Brezin et. al. 
and Wegner, and compare their density of states with various ways to reconstruct it 
from our computed moments. In section 6 we discuss the more general case of finite 
range impurity potential, still in the lowest Landau level. In section 7 we make a more 
restricted analysis for some higher Landau levels. 

2. Model and method of solution 



We choose units such that the magnetic length, — y^/e-B, becomes unity. I.e., we 
have (27r)~^ states per area in each Landau level. We measure energies relative to the 
(unperturbed) energy of the Landau level we project onto, and choose the energy scale 
such that the impurity potential becomes 



with g > 1. This is normalized to / ({fW{r) = 27r. The limit g ^ l'^ corresponds to a 
delta function potential. We assume a density //27r of Poisson distributed impurities, 
such that the random potential becomes 



G{x,y;t)^r{x,y)G{t), 





(1) 



V{r} = Y,W{f-Si), 



(2) 
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where the impurity positions {si} are independently and homogeneously distributed. 
I.e., there is on the average / impurities per state in a single Landau level. Since we 
neglect couplings to other Landau levels the Hamiltonian is simply equal to V, which 
we shall formally view as a perturbation. Thus, the zero'th order Green function equals 
the (kernel of the) projection operator onto the z/'th Landau level, 

GW(f,p;t)=n(f,p-). 

An explicit expression for this projection is 

where (x, y) and (^, rj) are the components of r and p respectively, and is the u'th 
Laguerre polynomial. For completeness we present a derivation of this well known result 
Appendix A| . Depending on the choice of gauge, the expression (^) may be multiplied 
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by a gauge factor exp [ir(r) — ir(p)]. Since is a projection it has the reproducing 
property, 

I dp P(fi,p)P(p,f2) = P(fi,f2). (4) 

Since L^{0) = 1, it follows from (^) that Vu{r,f) = {2tt)^^, in agreement with the fact 
that it must equal the number of states per area. 

A series expansion for the full G can be generated by repeated time integrations, 

fc=o ^■ 

We now average over the impurity positions {sj}. To this end we temporarily assume 
the system to be confined to finite area 2'kN with S = fN impurities, and let — > cxo 
afterwards. We find 



U{n,...,n) = V{n)---V{n) 



/ I[T.mn-s.J---W{n-sJ. (6) 

\i=l / j=l ij=l 

First consider the situation that all the Si^ correspond to different impurities (as 
is likely the be the case when all the positions fj are far apart). Then all integrations 
factorize, each giving a factor A^~^. Summing over all possible ways to pick k different 
impurities from a total of S, we get the contribution 

A^-'(^)=/(/-^)---(/-^) = /^'^-/' asAT^oo (7) 

to U{fi, . . . ,fk). The other extreme is to assume that all the Si- correspond to same 
impurity (as is likely to be the case when the density of impurities is very low, and all 
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Figure 1. The diagrams for the averaged Green function, up to 3rd order in 
the impurity potential. Open dots represent fixed coordinates, filled dots represent 
coordinates to be integrated over. A dashed line (with an implicit direction from left 
to right) represent the projection V^. A full line, or filled regions, represent various 
correlators of the random potential V. 
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Figure 2. The ith order diagrams for the averaged Green function. 



the fiS are close together). Integrating over the impurity position we get a contribution 
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to U{fi, . . . , fk). The full average must take into account all possible variations between 
the above two extremes. The complete expression is (as — oo) 

f/(fi, . . = E Ua,{f, . . . ,r ) ■ --U.^if,, . . . ,r ), (9) 

Q 

where the sum runs over all partitions Q of the set {ri, . . . , r/j} into j = l,...,k 
nonempty subsets of sizes ai,...,aj, and the f/^'s are symmetric functions of the 
coordinates in each subset. There are J2j=i S^^^ terms in the sum, where the S^^^^s are the 
Stirling numbers of the second kind[^]. Each term correspond to a Feynman diagram in 
the perturbation expansion. In our computer algebra program each fc'th order Feynman 
diagram is represented by a partition of the set {1, . . . ,k}. By generating all partitions 
of this set we obtain all A;'th order diagrams. They are all topologically distinct, and all 
occur with combinatorial factor unity. 

The 9 diagrams for the O'th through 3'rd order of expansion are shown in figure 0, 
and the 15 diagrams of 4'th order are shown in figure ^ Here open dots represent 
fixed coordinates, and filled dots represent coordinates to be integrated over. A dashed 
line (with an implicit direction from left to right) represent the projection Vu- The 
symmetric, translation invariant functions Uq{fi, . . . , fq) are represented by filled regions 
connected to q filled dots. For q = I they reduce to isolated filled dots, to each of which 
there is associated a factor /. 

According to (^^, a partition Q into j nonempty subsets gives a contribution 
which is proportional to . This is fully correct only in the limit A^ oo. A finite size 
correction is obtained by making the replacement f^^\ cf. equation 

The evaluation of the individual diagrams amounts to computing integrals like 

/ n 'Puir, pl)Vu{pl, P2) • ■ ■ V^{pk, p) Ua, (p., . . . , p.) ■ • ■ Ua^ (p., . . . , p.) 

1=1 

For the lowest Landau level, = 0, this is just a sum of Gaussian integrals. In the 
higher levels it becomes a sum of Gaussians multiplied by polynomials. 

Since the averaging process restores translation invariance, each diagram must be 
a translation invariant (with respect to the magnetic translation group) expression 
constructed only out of states in the z/'th Landau level. This forces it to be proportional 
to Vu{r,p). To show this more explicitly, we start with the most general expansion of 
an averaged diagram, 

Gair, P)= j dpdg C(p, q) ipp{r) ^g(p)*. 
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where {ipp} form an orthonormal basis for the z/'th Landau level. 

In the Landau gauge, = 0, we may choose this basis such that the magnetic 
translation group acts as 

t{ax) : Vfc e**""' ipk, t{hy) : i^k+b', 
where a' is proportional to a, and b' is proportional to b. The requirement that 



m)GQ] {f,p) = GQ{f,p) 

implies that C{p, q) = C{p — q), and the requirement that 



[t{ax)GQ] {f,p) = GQ{f,p) 



implies that G{p — q) = Go6{p — q), with Co a constant. Thus, GQ{f,p) oc Vy{r,f>). 
A graphical representation of this relation is illustrated in figure |] (except that the 
split-off "vacuum diagram" should not be counted with its usual combinatorial factor). 
As a control of our computer algebra program, we have verified this relation explicitly 
for a large number of cases. Using it reduces the evaluation of /c'th order diagrams to 
calculating determinants oi k x k matrices (when employing complex coordinates in the 
lowest Landau level). 

Comparing the expansion 



-t) 



k 



G(f, p; t) = P,(f, p) Y: '-jr Gk = V^if, p) Git) 

k=0 

with the general expansion in energy eigenfunctions 
G(f, p;t) =J2 ^c>{r) C(P) exp i-tEa) , 

a 

we obtain after setting p = f and integrating 

N G{t) = exp i-tEa) = N f dE D{E) exp {-tE) , (10) 

where D{E) is the averaged density of states (normalized to / dED{E) = 1. I.e., G{t) 
is simply the Laplace transform of D{E), and the information obtained from the fc'th 
perturbative order is exactly the fc'th moment of the averaged density of states. 



Gk = jdED{E)E\ (11) 

3. Counting and resumming graphs 

We define 

G{u) ^ - / dt e*/" G{t) = / dE D{E) — = YGk u'+\ (12) 

Jo J 1 — Eu ^ 
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Figure 3. Due to translation invariance, every diagram contributing to the averaged 
Green function becomes proportional to the projection onto the given Landau level. 

where the series should be interpreted in an asymptotic sense. The integral converges 
for u < 0, and G{u) can be extended by analytic continuation. As noted, the graphical 
expansion illustrated by figures amounts to expanding G into a formal power series 
in u. There are various ways to resum graphical expansions of this type. The simplest 
and most useful one is to sum up sequences of V^s connected by isolated dots (i.e., 
factors of Ui). This amounts to making the replacement jzj^j, and summing 

over graphs containing no isolated dots. The latter sum defines a function G{u), related 
to G by 

G{u) = G{—^). (13) 

The fc'th order graphs for G{u) ~ Y^'kLoGku'^^^ are in 1-1 correspondence with the 
partitions of the set {l,...,k} into subsets of sizes aj > 2. The number of such 
partitions is significantly smaller than the total number. For instance, inspection of 
figures |l]-0 reveal that only the graphs (0), (1), (26), (3e), and (4£)-(4o) contribute to G 
up to the 4'th order. The total numbers of graphs which contribute to G{u) {Gk) and 
G{u) {Gk) are listed in table |I|. 
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Table 1. The number of graphs contributing to various functions. Gk- The total 
number of graphs contributing to the Green function. Gk'- The number of graphs 
not containing isolated dots. Tik'. The number of 1-particle irreducible graphs not 
containing isolated dots. Sk- The number of skeleton graphs. Ck- The number of 
irreducible terms in the answer. 



k 


Gk 


Gk 




Sk 


Ck 


1 


1 








1 


1 


2 


2 


1 


1 


1 


1 


3 


5 


1 


1 


1 


1 


4 


15 


4 


3 


2 


2 


5 


52 


11 


9 


6 


2 


6 


203 


41 


33 


21 


7 


7 


877 


162 


135 


85 


10 


8 


4140 


715 


609 


385 


36 


9 


21147 


3 425 


2 985 


1907 


85 


10 


115 975 


17 722 


15 747 


10 205 


319 


11 


678 570 


98 253 


88 761 


58 455 


1 113 


12 


4 213 597 


580 317 


531561 


355 884 


5 088 


13 


27 644 437 


3 633 280 


3 366 567 


2 290 536 




14 


190 899 322 


24011 157 


22 462 017 


15518391 




15 


1 382 958 545 


166 888165 


157 363 329 


110 283179 




16 


10 480142147 


1216 070 389 


1 154 257 683 


819 675 482 





Another standard resummation is to define through the relation 

G{u) = (14) 

where only 1-particle irreducible graphs contribute to S. We may combine this with the 
resummation of isolated dots above, to obtain 

S(«) = /M+(l-/n)S(- -). (15) 

1 - /m 

Here S is the sum of all 1-particle irreducible graphs which do not contain isolated dots. 
The numbers of such graphs are also listed in table |1|. As can be seen, there is some 
reduction in the number of graphs to be calculated, but not a very significant one. Of 
the graphs in figures the only additional saving is the elimination of graph (4£). 

A further reduction in the number of contributing graphs is obtained by first 
summing the up the skeleton diagrams S{u). This is the set of graphs such that their 
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corresponding "vacuum diagrams" (cf. figure |^) are 1-particle irreducible. This function 
is related to S by 

7/ 

S(..) = [l-S(t.)]5(^-^). (16) 

(The function S{u) = S{u) — fu is related to S('u) in the same way.) The numbers of 
graphs contributing to S are also listed in table |^. Expanding S to second order in its 
argument, and solving the resulting 2nd degree algebraic equation for E (or, equivalently, 
G), constitutes to the so-called self consistent Born approximation (SCBA). Of the 
graphs in figures |l]-0, the only additional saving when restricting to skeleton diagrams 
is the elimination of graph (4m). 

One way to employ the resummation methods above is to generate a smaller set 
of diagrams (e.g., for G{u), £(«), or S{u)) to a given order, and then find the full 
perturbation expansion to the same order by use of the algebraic relations above. 
However, the benefits of going beyond G{u) are marginal, since the reduction in number 
of graphs is rather small while more computer time is needed to classify graphs. We 
have in our computation evaluated all graphs contributing to G{u). As a check on the 
computer algebra we have in addition summed the subsets of these graphs contributing 
to S(m) and S{u), and verified that they reproduce the same end result for G{u). The 
simplest check of the program is to verify that it actually generates the number of 
graphs listed in table |l], since these numbers are computed in an entirely independent 
way (starting from the Stirling numbers). 

Another way to use the resummation methods above is to combine a finite 
perturbation expansion from S(n) or S{u) with the exact algebraic equations, thereby 
generating some approximate, but infinite, series for G{u). Considering the numbers in 
table ID, we note that the resummation of isolated dots captures a large percentage of 
the total number of graphs. However, the physical effect of this resummation is rather 
uninteresting — it merely corresponds to an overall shift in the energies. Of the remaining 
graphs we note that the majority of them are complicated ones (i.e., skeleton graphs) 
which cannot be generated by resumming lower order terms. Thus, the resummations 
above underestimate the fc'th perturbative order by increasingly large amounts as k 
increases. In the absence of any physical reasons for selecting a particular class of 
diagrams, we believe it is much more sensible to extrapolate the perturbation series in 
a statistical sense, by considering how the total number of graphs, and their average 
value, varies with the perturbative order k. 

4. Calculated moments 

We have calculated the perturbation series up to the 12 'th perturbative order, for general 
values of impurity density / and potential range g. Each graph contributes a rational 
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function in g, multiplied by some power of /. The general answer rapidly becomes too 
complicated to present here. We present the result in terms of the moments Gk of the 
density of states. The general expressions for the first five moments are listed in table 0. 
The full result for the case of 5-function impurities, = 1, is listed in table |^. We list 
results for some additional values of g (in floating point form) in [Appendix B. 



Table 2. The first five energy moments for general g, in the lowest Landau level. 
k Gk 

1 / 

2 ^J + f 

A 2 f I 25 g*+25 g^+2 /•2 , 6 j-3 , f 4 

t- 16 f I 210e"'+420 e'' + 170 r2 I 85 g;*+85 g'^ + W f 3 , 10 f 4 , fb 

T^TWpTi^ ^ e(3e^+i)(3e^+5)(e^+i)-' ^ e^(e^+i)(3 g^+i)J g ■' ^ 

U 2*° r (k\ (9fc-ll)(e'^ + g") + (2fc-6) rk-2 , fc(fc-l) jfc-l , ffc 

(e+i)''-(e-i)^*' ^ V3y 4e^(e^+i)(3e^+i) 2g J ^ 



It is an amusing exercise to search for general patterns in table ^ and table ^. The 
coefficients of the highest powers of / are straightforward to find, since the series for 
G{u) up to the 2j'th perturbative order determines all the terms proportional to f^, 
f^~^, . . . , /^~"' in the /c'th perturbative order. The lowest powers of / are of more 
physical interest. Consider first the case of 5- function impurities. The general pattern 
of the /^-terms looks a bit complicated, but it is easy to verify that it fits the formula 

Lfe/2J 

0(f). (17) 



Lfe/2J , / 

G, = f + f E - 




These moments are reproduced to order by the density of states 

DiE) = (1 - /) 6iE) + \ f e{E) 9{2 -E)\E- (18) 
The physics behind this expression is as follows: 

(i) Consider S = f'N < N arbitrary placed impurities in a single Landau level, for 
which the Hilbert space of wave functions is A^-dimensional. The condition that 
a wave function vanishes on all impurities imposes 5* constraints, for which the 
solution space is at least {N — S") -dimensional. These solutions will have exactly 
zero energy, which is also the lowest possible energy. Hence, there must for all 
/ < 1 be a (5-function contribution to the density of states, of strength (1 — /). 



Table 3. The twelve first energy moments for (5-function impurities {g — 1). 
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+ f 
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+ 3P + 


f 




4 
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+ f/^4 


6/3 + /4 




5 


/ 


+ f/^4 


fP + iop + f 




6 


/ 


1 137 f2 


f 2^/3 + 115 /4 + 15/5 + /6 




7 


/ 


1 245 r2 


f2fiI/3 + lf'/4+ 245^5 + 21/*^ + /^ 




8 


/ 


1 871 r2 

+ TTJ ' 


f + ^T-P + 749/5 + 231/6 + 2%f + p 




9 


/ 




f 47443 ^3 + 29091^4 ^ 31121^5 + 1348^6 + 399^7 + 35^8 


+ P 


10 


/ 

+ 


1 4629 f2 
^ 20 

45/9 + / 


1 241715 j-3 1 452455 f4 , 292901 j-5 , 98261 f6 1 8085 j-7 , 
'72''' 36-'' 16-'' 8-'~'~2-'"'" 

10 


645/8 


11 


/ 


1 8349 f2 
20 


1 3068197 /?3 1 3047209 j-4 , 11665093 f5 , 1163173 fQ , 267377 r7 
360 ' 72 ' 144 J ' 16 •' ' 8 




+ 


16137^8^ 


-990/9 + 55/1" + /" 




12 


/ 
+ 


1 45517 J- 
60 J 
1940433 f7 
8 J 


2 1 5201203 f3 , 303556067 f4 , 16554385 f5 , 114943133 f6 
240 ' 2160 J ' 48 ' 288 

+ 649^231^8 _^ 149 82/9 + 2915/1" + 66/" + /^^ 





At very low / the impurities are very far apart. Each of them will bind one state 
of energy E = 1. Hence, they will contribute a term f 6{E — 1) to the density 
of states. This is sufficient to reproduce all moments to order /. To order p we 
must consider the mixing between states at different impurities (to this order only 
mixing between two impurities). The energies for two impurities at distance r are 
E± = 1 ± e"*^^/^. Integrating over the distribution of relative distances we get a 
contribution to the density of states, 



AD{E) oc y df [6{E - 1 - e""'/^) + - 1 + e" 

/ — [6{E~l-x)+6{E-l + x)] 
Jo X 



oc 

oc e{E)e{2-E) \E-1\~\ 

This distribution is not normalizable, but has finite moments {{E — l)'^). The lack 
of normalization is due to the fact that additional impurities must be taken into 
account at very large separations r (of order /~^/^). A simple way to model this is 
by introducing a distribution dp{r) = e'^^"^^ kfdr'^ for the relative pair distance (the 
total density of pairs being \ f)- The requirement that all moments be reproduced 
to order p fixes k = \. This in turn leads to (|18]). 
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(iii) If one repeat the previous argument with the finite size correction, f^"^^ = 

f{f — N^^), one is lead to the conclusion that the exponent in (18) should also 
undergo a finite size correction, 

\E - ^ \E - if . 

It is also easy to understand the physical origin of the low-/ behaviour of the general 
moments, 

Gk = , ^^.f , T^/ + 0(/'). (19) 

{g+ l)" - {g- l)" 

This is related to the behaviour of electrons near the single impurity potential (|l]). With 
the impurity at the origin we find eigenstates (in the symmetric gauge, and cylinder 
coordinates) 

V>,(f) = (27r2^£!)-i/2 /g-iVe-'-V4^ 
and corresponding energies (cf. Appendix A) 

^-(jTrjlfri)'- (20) 

With a small fraction / of such impurities the contribution to the density of states 
becomes 

D{E) = fJ26{E-E,). (21) 

£ 

The sum must be cut off when the extension of the wave function becomes of order 
/~^/^, i.e. for i of order /^^. This is important for obtaining a normalized D, but can 
be ignored when calculating its moments. Thus we get to order / 



oo / 2 

G, = fY: = f [—TT 



\g+l 



k 



-1 



4-1. The cumulant expansion 

The information in table ^ and table § can be compressed somewhat by rewriting it 
in terms of cumulants. The moments (E^) are related to the Laplace transform of the 
density of states, 

oo f ^•jfc 



G(t) = y^ di?D(E)e-*^ = 5:^G,. (22) 

By rewriting 



k=0 



Git) = exp (x(t)) = exp i^Y. Xkj (23) 
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Table 4. The first five cumulants for general g, in the lowest Landau level. 
k {-!)'' klxk 

1 / 

2 ^/ 

A 2 r _ 1 f 2 

c 16 f 80 f 2 

5e*+i0e^+i-' (3e^+5)e(3 g'^+i)J 



2- 

(e+i)'=-(e-i)'=-' 



Table 5. The 12 first cumulants for 5-function impurities, in the lowest Landau level. 
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-if 


8 


f- 


653 j-2 1 
12 J 


757 j-3 109 j-4 

6 8 


9 


f- 


^/2 4 


11603 J-3 2085 j-4 
18 8 


10 


f- 


5591 J-2 
20 


1 204725 J-3 17065 j-4 , 2971 /-5 
^ 72 6 ^ 16 ■' 


11 


f- 


12111 r2 
20 J 


1 820787 ^3 841489 j-4 , 87197 f5 
^ 72 J 36 ' 16 


12 


f- 


77303 J-2 
60 ■> 


1 10270399 J-3 69677951 f4 , 87197 j-5 124513 j-6 
^ 240 •' 432 ■' ' 16 •' 32 ■' 



we obtain the cumulant expansion. No information is lost when going between the set 
of cumulants {Xk\k = 1, . . . , m} and the set of moments {Gk\k = 1, . . . , m}. 

We show the five first cumulants for general values of / and g in table |, and the 
twelve first cumulants for (5-function potentials {g = 1) in table ^. Note that a factor 
(— l)'^//c! has been split off in both tables. The cumulant expansion may also be viewed 
as yet another way of resumming a perturbation expansion. We return to this below. 
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5. Comparison with the exact results of Brezin et. al. and Wegner 

For the case of 5-function impurities, ^ = 1, exact formulas for the density of states have 
been found by Brezin et. a/.0], and by Wegner P]. The result by Wegner correspond to 
taking the additional limit of / — > oo. Their expressions for the density of states are 
highly non-trivial, and they provide very useful checks on our results and methodology. 

5.1. The result of Brezin et. al. 




Figure 4. A 3I?-plot of the density of states found by Brezin et. al., as function 
of energy E and impurity density /. In addition to the plotted density there is a 
contribution (1 — /) 5{E) when / < 1. 



The results found in reference |0] specializes for our model (with q = 1) io the 
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integral 

Db(E) = - Im— In / dt expliEt - f I{t) ] , where 

TT OE Jo 

m = jl^{l-e-^)=it + \e-i,t^-.... (24) 

It is a rather challenging task to evaluate this expression numerically. As was pointed 
out in [0, Db{E) has a lot of interesting structure, in particular for low density of 
impurities /: 

(i) For / < 1 there is a 5-function contribution to Db{E), 

DBiE) = il-f)6{E) + ---. (25) 

As explained before (cf. the discussion around equation ([l8|) ), this is due to the 
fact that we may arrange for a fraction 1 — / of all wavefunctions to vanish on all 
impurities. This phenomenon was already pointed out by Ando^]. 

(ii) As E ^ 0^ one finds the behaviour 

I E~f < / < 1, 

DBiE)r^\ E-^/\og\E) / = 1, (26) 
i E^-' 1 < /. 

The increasingly singular behaviour as / increases from to 1 is due to zero-energy 
states moving out into the low E region. As for the large /-dependence, a crude 
qualitative understanding is obtained by considering an arbitrary placed, maximally 
locahzed state. 

The probability Prob(ii^ < Eq) that this state has an energy less that Eq is equal 
to the probability that the distance r from r'o to its nearest impurity satisfies 
U exp(-rV2) < Eq. Now, 

Prob (exp(-rV2) < Eq/2txU) = (^o/SttAT)^ , 

which predicts D{E) ~ E^~^. This is an underestimate (because our assumption 
that the state was maximally localized, and arbitrary placed, cannot be expected 
to hold), but the dependence on / turns out to be correct. 

(iii) As ^ 1 one finds a diverging density of states 

Db{E) r^\E- for < / < 1, (27) 
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and this singularity continues as a cusp singularity for 1 < / < 2. As discussed 
before (cf. equation ([18|)), this singularity originates in states localized on top of 
single impurities. The interaction with nearby impurities leads to a level broadening 
around E = 1. 

(iv) As E ^ 2 there is a cusp singularity, 

^Db{E) ^ \E - 2\f~\ forO</<l. (28) 

This singularity is due to the fact that a state must interact with at least three 
impurities to have an energy E > 2, while two impurities are sufficient for E < 2. 
By a similar reasoning one can also understand why there are (increasingly weak) 
singularities in D^{E) at = 3, 4, . . .. 

(v) Db{E) is a fairly smooth function for / > 2. 




Plots for the integrated density of states, dE'DsiE'), where given in g] and 0. 
The plots of Db{E) itself look more interesting; we give examples in figures We 
shall return to the a reconstruction of Db{E) from our calculated moments, but first 
consider the limiting case of / — oo. 
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5.2. The limit of high impurity density, and the result of Wegner 

Consider now equation (24) for very large /. The main contribution to the integral will 
come from small t, hence we may approximate 

i{t) ^it + \e. (29) 

This leads to the expression 

fD^{E)^D^{t) = -Im|-ln(l + erf(ie)) 

^ J_ exp(6^) ^ J_ exp(e^) 

vr3/2 1 _ erf2(ie) vr^/a i + irf ^(e) ' ^ ' 

where e = {E — f)/ a//, and irf(e) = — i erf(ie) is a real function of a real argument. This 
is equivalent to the result by Wegner f^. 

The Fourier transform of -Dw(e) has the series expansion 

/oo 
dee-^'^^Dw(e) 
-oo 

= l-l(e^)c.^ + i,(6^)c.^-i(e^)..^ + i(e«)..«-.--, (31) 
with (e") = de Dw(e) e"- One may verify (numerically) that 

(e^) = 1, (6^) = I, (e«) = f , = f , (e^°) = ^, (e^^) = (32) 

All odd moments vanish because -D-w(e) is an even function of e. 

Now turn to the problem of determining D{E) from e.g. the cumulants in table ^. 
Formally this may be done by inverting the Laplace transform ( p2l - p^ ). With t = iu this 
becomes 

1 /.oo , ^ °° 

D{E) = — / dwe^(^-^)^ expJ2im)'xk. (33) 

ZTT J-oo f,^^ 

To investigate the large-/ limit we introduce e = {E — f^j^f] ■, ^ = "^vT; cind keep only 
terms which remain when / ^ oo. Thus, in this limit, 

fD{E)^— r duj e'^'^ D{uj), (34) 



27r J- 

where we find from table ^ 

f)(, ,^ _ pYr. r 1 , ,2 1 4 7 , ,6 109 , ,8 2971 , ,10 124513 , ,12 \ (o^\ 

One may verify that Dw('^) = -D(u;) to the order in u we have computed. 



5. 3. From cumulant expansion to density of states 

We now construct approximants Dn{e) to D{e), using equation (p3D with the cumulants 
X2, ■ ■ ■ ,Xn included in the sum. The result is surprising and instructive. To investigate 
convergence we evaluate the densities at e = for increasing n. We find 

Dn{Q) = (0.3788, 0.3744, 03721, 0.3710, 0.3703) for n = (4, 6, 8, 10, 12). 



19 

This sequence converges nicely with n (see figure fitting the formula 

D„(0) ^ 0.3675 + 0.02456 n'^ + 0.09745 + ■ ■ ■ . 

However, when extrapolated to n = oo we obtain Z}oo(0) ~ 0.3675, which is about 2.3 
percent higher than the exact result, 

Dw(0) = 2/7r3/2 = 0.35917.... 



0.38 



o 

e 0.37 



0.36 



0.1 , , 0.2 0.3 

1/n 

Figure 6. The approximants Dn(0) as function of together with a fitting function. 
It is clear that the approximants converges to a value which is different from the 
(indicated) exact result, i3w(0) = 0.35917. . . 

The situation is further illustrated in figure |^, where we compare the functions 
Dio(e), -Di2(e) and Dy/{e). The curves for -Dio(e) and -Di2(e) are almost 
indistinguishable. This indicates that the sequence of approximants has converged. 
However, the limit function turns out to be negative for certain ranges of e. Thus, even 
without knowledge of the exact result we could have concluded that something was 
wrong. 

The origin of the problem can be understood by investigating the behaviour of 
Dt^{uj). This is plotted in figure ^, together with the approximants D^^uj), Ds{uj) and 
Di2{uj). We observe that Dy^^uj) becomes negative for some Ic^l > cuq ~ 2.3. Since all 
the Xk are real, the cumulant expansion can never reproduce such a behaviour. In fact. 
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0.4 












1 



2 



3 



4 



e 



Figure 7. A comparison of the approximants Dio(e) and -Di2(e) (dashed lines) with 
the exact result _Dw(e) (fully drawn line). Since all curves are symmetric with respect 
to e ^ — e only positive e is plotted. Dio and D12 are almost indistinguishable, 
indicating good convergence. However, the approximants have converged to a result 
which is different from the exact one. 



all terms in the exponent in equation (^) seem to have the same (negative) sign. Thus, 

we expect the series in the exponent to converge to \ogD^J\/{u) for \u\ < Uq, and diverge 
to —00 for |u;| > uq. Thus, our approximants Dn{u!) will converge to the (wrong) limit 



The curves in figure |^ confirm this behaviour. The lessons of this section are that (i) 
arbitrary resummations may lead to misleading results, and (ii) even if a sequence of 
approximations converges, it may converge towards the wrong result. 

5.4- From moment expansion to density of states 

In the previous section we realized that the cumulant expansion for D^^{uj) has a finite 
radius of convergence, c^o ~ 2.3. Here we shall first show that the moment expansion 
for Dt^{uj), cf. equation (31), has an infinite radius of convergence. To this end we 
must evaluate (e") as n — ^ 00. This quantity will receive its main contribution from 
the region of large |e|, where we may approximate -Dw(e) ~ 27r^^/^ exp(— e^). An 





for \uj\ < ujq, 
for \uj\ > Uq. 



(36) 
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Figure 8. Fourier transforms of density of states. Tlie exact result Dw{lo) (fully 
drawn line), is compared with the approximants r>4(a;), -Dg(a;) and £'12 (w) constructed 
from the cumulant expansion. The exact function becomes negative in some regions, 
a property which cannot be reproduced by the approximants. 



asymptotic evaluation of the integral de e"+^e then reveals that it receives its 
main contribution from |e| ~ ^1 + n/2 (confirming the approximation), and that 

It follows that the series D^{uj) — Y^^=q{—i-^'^Y (e^")/2n! converges for all uj. Clearly, 
a term-by-term integration of this series will not lead to a sensible density of states. 
However, wc may extract a convergence factor e~"'^^/^, and series expand the quantity 
e"'^^/^ Z)vv(^)i which also has an infinite radius of convergence. With (the somewhat 
natural choice of) a = 1 this method gives better results than the cumulant expansion, 
but the convergence with n is unimpressive. (With a priori knowledge of the answer it 
is indeed possible to choose the parameter a such that an excellent reconstruction of the 
density of states is obtained. However, we have not found a good objective criterion for 
choosing a, which would work in more general circumstances.) 
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5. 5. The known moments as constraints on the density of states 




Figure 9. A comparison of the exact density of states, Db{E) (fully drawn lines) 
and the density D{E) reconstructed from the 12 first moments (dashed lines) for a) 
/ = 1.99, b) / = 2.5, c) / = 3.0, d) / = 3.5, e) / = 4.0, g) / = 4.5, and h) / = 5.0. 



Given only a finite set of moments, {E'') ior k — 0, . . . ,n, any non-negative function 
D{E) which reproduces them is in principle a possible solution for the corresponding 
density of states. In practice, D{E) may be expected to have some simple asymptotic 
behaviour for large E. Since {E^) for large k is mainly determined by D{E) for large 
E, we may use the asymptotics of {E^) to estimate this tail of D[E). Further, for 
the models at hand, we know that D{E) = for £J < 0, and we may have some 
independent information about how D{E) behaves at small E. With the behaviour of 
D{E) constrained from both sides, one expect it to be well determined by its lowest 
moments — provided it is a reasonably smooth function. 

This strategy gives a satisfactory method to reconstruct the density of states. We 
find that the asymptotics of our calculated moments fits well to the behaviour 

Such a behaviour is reproduced by a density of states which behaves like 
D{E) ~ E" exp (-aE'^ + hE) as E ^ oo. 
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Figure 10. The accuracy of the reconstructed density. We plot the quantity 

{D{E) - Dji{E)) / d'^^""''^ for / = 3 (fully drawn Hne) and / = 4 (dashed line), where 
jjimax) ^j^^ maximal value of the function Db{E). 



We have the connection a = (2a)^^/^ and P = b{2a)^^. The index c does not enter 
to the order we consider here. With the coefficients a and b determined we write the 
density of states in the form 

D{E) = Ef~^ exp (-aE"^ + hE) x P„(E), (37) 

where P„ is a n'th order polynomial We determine its n + 1 coefficients such that all 
the moments {E^), . . . {E^) are reproduced. This construction leads to a D{E) which 
compares favourably with the corresponding Db{E), at least when / is somewhat larger 
than 2, cf. figure ^ 

Since the difference between the exact and the reconstructed curves is mostly 
invisible in figure ^, we show some examples of this difference in figure |10|. The 
accuracy is seen to improve with increasing /, most likely because the function we try to 
reconstruct becomes smoother (thus this trend may not continue to arbitrary high /). 
In this reconstruction we have built in the known low-E behaviour, D{E) ~ E^'"^. We 
have noted that the overall reconstruction is fairly insensitive to shifting the exponent 
away from / — 2. Choosing the correct exponent may give somewhat better convergence 
with n. 
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6. Density of states for general q 

The method used in the previous subsection can be used equally well to reconstruct the 
density of states for q > 1. The main difference is that we have no a priori knowledge 
of the \ow-E behaviour of D{E). However, in the previous section we gave a heuristic 
probability argument for the /-dependence of the exponent. This argument may be 
repeated for a potential of finite range. It suggests that one should make the replacement 



f^f=l{l + g)f (38) 

in the exponent when ^ > 1. Thus, we have imposed the requirement that D{E) ~ E-l'~'^ 
as — 0^. The resulting density of states is shown in figure for a set of f)- values. 
We now longer have exact results to compare against, but comparison with numerical 
simulations show excellent agreement. We believe the difference from the exact curves 
would not be visible in the plots. 




Figure 11. Constructed density of states for / = 4 and a) g = 1.0, h) g = 2.5, and c) 
£> = 5.0. 



As g increases with / fixed, the distribution approaches a Gaussian of width y f / g, 
centred around E = f: 

D{E) ^ {27rf/g)-'/'exp ("^^^^7^) • (39) 
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In terms of the variable e = {E — f)y g/ f, the first high-f) correction to this distribution 
is determined by the third moment 

(e^) ^ liQf)-'/'. (40) 

In general, the fc'th order correction (as e.g. defined by the fc'th cumulant in e) vanishes 
like p(2-fe)/2 as — >• oo. 

7. Higher Landau levels 

Calculation of averaged Green functions in the higher Landau levels proceed as in the 
lowest one. The changes are that we measure energy relative to the band centre, 
E = (z/ + 1/2), and use the appropriate expression for the projection operator, i.e. 
V ^ Vu, where u is the Landau level index. The (integration kernel) of the projection 
operator for arbitrary Landau level is given in [Appendix A| . The perturbation expansion 



for the averaged Green function continues to be equivalent to computing the moments 
of the energy density of states. Since we shall only perform a low order calculation "by 
hand" we may as well calculate the moments directly. (The main advantage of working 
with the Green function instead of directly with the moments is that we don't have to 
worry about combinatorial factors. This reduces the time it takes to code, debug and 
run computer algebra routines.) 

7.1. Evaluation of the three lowest moments 



The diagrams for the first three energy moments are shown in figure ^ . 
We have 

{E), = f, {E')u = f, {E')ui = f, {E'^)i2 = 3f{E')2, (41) 

which reduces the problem to evaluating {E'^)2 and {E^)i2. These quantities depend on 
the Landau level index, which we shall indicate in our notation below. The evaluation 
of {E^)2 reduces to computing 

where is the z/'th Laguerre polynomial. Since this is a Gaussian times a polynomial, 
it is simple to evaluate for the first few u. We find 



2\H 



2 



-J 

1 

2 (g-l)4+4(g_l)2 + l 

g5 J 

3 (g-l)«+9(g-l)^+9(g-l)^+l r 
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{E'h 



{E'y. 



{E%i 




Figure 12. Graphs constituting the first three energy moments (E), {E"^) and {E^). 
Note that the statistical weight is a factor fc! higher than the usual combinatorial factor, 
where k is the order of the moment considered. 



This fits with a general expression 



n2\(i') _ / ^-(2i/+l) ^ 
fc=0 



(42) 



where is a rescaled, effective form of g, allowing us to write the second order moment 
for arbitrary u in the same form as the second order moment for the lowest Landau level. 



We show in figure ITB| how this quantity varies with p and u. 
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A similar evaluation of {E^)y for the first few u gives 



1 (iTW (^^^^ ~ -^^^^^ + ~ ^^^^^ + ^^^^ + 84^ - 11) / 

2 (x+l^ (729^12 _ 5832^11 + 24786^^° - 61560^^ ^ 94527^8 _ 72144/ 

-23076/ + 114192/ - 104697/ + 29976/ + 11826/ - 4632^ + 1) / 



Here the generalization to arbitrary u is not apparent, except that specializing to ^ = 1 
gives (-E'^)3'^^ = /. The expressions above are related to the cumulants by 

X^:^ = {E^^:^ for n = 1, 2, 3. (43) 

77". 

(This relation does not generahze to higher n.) 
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1.2. From moments to density of states 




E 

Figure 14. Density of states constructed from {E), {E'^), and (E^) (dashed lines) for 
a) v = 0, h) v = I, and c) v — 2, with f = g = 2. For comparison we show the density 
constructed from the 12 first moments in the lowest Landau level (full line). 



In this case we have calculated too few moments to use the method from sections 5-6. 
Instead we parametrize the density of states as 

D{E) = CE' exp (-aE^) (1 + bE) , (44) 

where C is a normalization constant, and the parameters a, b, c are chosen such that 
{E), {E^), and {E^) are reproduced. To get some measure of the accuracy obtained, 
we also use the same procedure for the lowest Landau level. For f = g = 2 we find the 
following fitting parameters 



u 


a 


b 


c 





0.22174 


0.06732 


1.10221 


1 


0.47681 


0.22997 


2.98668 


2 


0.64008 


2.50640 


3.77993 



The resulting densities are plotted in figure 14. For comparison we also show the density 
found for the lowest Landau level when all 12 moments are used. It seems that D{E) 
can be reproduced fairly well from only three moments. 
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8. Concluding Remarks. 

We have in this paper generated and analysed a high (12'th) order perturbation 
expansion of the averaged Green function, in a class of models inspired by the integer 
Quantum Hall Effect. The generation process can be fully automatized using computer 
algebra programs. The order which can be obtained by our programs is mostly limited 
by available CPU time, as the number of diagrams exhibit factorial growth with the 
perturbativc order. With present technology about 3 more orders would be feasible for 
the quantity considered. 

For the models considered the n'th order perturbation expansion provides exact 
information of the n first moments of the density of states versus energy, D{E). When 
D[E) is a reasonably smooth function it can be reconstructed to high accuracy (better 
than 1 percent) from the computed moments. The results are also useful to check and 
complement other approaches, like numerical simulations on finite size systems. 

There are several directions to extend our method. The one of most immediate 
interest for the integer Quantum Hall Effect is to extend the analysis to the averaged 
2-particle Green function. This quantity encode information about transport properties. 
We have done some initial investigations in this direction. It is fairly straightforward 
to implement a computer algebra procedure which automatically generate all diagrams, 
and evaluates them. The number of diagrams will be about one order of magnitude 
larger (for e.g. the 12'th perturbativc order). The information generated can be viewed 
as various moments of a 3-variable spectral function. It is not clear how easy it will be 
to reconstruct this spectral function, or the interesting physical information it contains, 
from the information generated by the perturbation expansion. 
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Appendix A. Projection operator onto the r/'th Landau level 

With appropriate units and gauge choice, the imaginary time propagator for a 2D 
electron (with charge — e, where e > 0) in a magnetic field solves the initial value 
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problem! 



d 



G{f, p; t) 



dr 
G{f,p;0) 



' d i V fd i V 



G{r, p; t) 



5{f-p) 



(Al) 
(A2) 



By making the ansatz that G is a Gaussian in (x, ^, rj), with prefactor and coefficients 
which depend on r, one obtains the solution 
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1-e 



(A3) 



where e = e and R = |(r — p)^. Since the energy of the z/'th Landau level is + ^ 
u = 0, 1, . . ., an alternative expression for G is 



G(f,p;r)=e--/2^ V,{f,p)e- 

u=0 



(A4) 



where is the projection onto the i^'th Landau level. Thus, by series expanding 
equation (|A3|) in powers of e we easily find the (integral kernels of the) projection 
operators onto the various Landau levels. The ffist few examples are 



Viif,p) 



2tt 
1 

2^ 



R/2 



V2{f, p) = ^ e'^-v-iy)/2 (1 _ 2i? + e-«/2. 
zvr 



The general expression is 



1 



V,{r,p) = ei(-^-«^)/2 L,(i^)e-^/^ 
/vr 

with L,y the z/'th Laguerre polynomial. 



(A5) 
(A6) 
(A7) 

(A8) 



Appendix A.l. Eigenvalues in a rotation symmetric potential 

The projection operators can be expanded in a basis of rotation eigenfunctions. With 

X = r cos if,y = r sin (f,C,=r' cos (f', rj = r' sin (f': 



Pu{r,p)= E e-''^^~^'^ Ft\r) Ft\r') 



(A9) 



f Here r and p are 2£) vectors with components (cc, y) and (^, 77) respectively. Units for length and 
time are chosen such that £b = {h/eBY^'^ and ujb = eB/m equals unity. Orientation is chosen such 
that B points along the positive z-axis, and the problem is formulated in the (radial) symmetric gauge, 
f • A{r) = 0. 
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If we introduce a rotation symmetric potential, W{r) = W{r), {e ^^'^ Fi{r) } continues 
to be a basis of energy eigenf unctions. The corresponding eigenvalues are 



= 2n r rdrF^''\ryW{r). 
Jo 

Combining ([A9[ - |AlCl| ) gives a generating function for the eigenvalues, 

^0 „ 



(AlO) 



(All) 



e=~u 



Here p is chosen such that r' = r and ip' = 0. With W{r) as defined in equation (|I|) the 
integration is simple, and we get a generating function 



^ + 1 



{ae'^^y (1 - be-'^^y (1 - ae 



(A12) 



where a = {q — \) / {q + 1) and h = {g — ?>) j {q — 1). It is easy to verify that 

oo 

ZM(1) = £ Ef^ = 1. 
Series expansion gives the eigenvalues 



J+2p 



g + l 



E 

k=0 



kj [u + i-kj [a 



(A13) 



We find that ^ ^ 6e, SuS g 1+. More explicit expressions for the first three Landau 
levels are 



(0) 



2 ( g-V 

g + l U + 1 



E, 



(1) 



(2) 




4(£ + l) + (f?-l)^ 



(A14) 
(A15) 



^ + 1)5 \^^+l 

8{i + 2){i+l) + 8{i + 2){g- 1)2 + {g- 1)4 



(A16) 



Appendix B. Tables of cumulants. 

We here present the 12 first cumulants for some selected values of g. The exact 
expressions are much too long to present in full (except when g = 1), thus we give 
only their numerical approximations. It is apparent from these tables that the higher 
cumulants vanishes quite rapidly when g is large (and / is of order one). 



Table Al. The 12 first cumulants for g = 1. 



k 


(- 








1 


/ 








2 


/ 








3 


/ 








4 


/ 


- 0.500000/2 






5 


/ 


- 2.500000/2 






D 


7 


— O.iDDDD// 


_i_ 1 ^Kr\r\r\r\ f 3 

-\- 1. i OUUUU/ 




7 


/ 


- 22.16667/2 


+ 19.25000/3 




8 


/ 


- 54.41667/2 


+ 126.1667/3 


- 13.625000/* 


9 


/ 


- 125.7500/2 


+ 644.6111/3 


- 260.62500/* 


10 


/ 


- 279.5500/2 


+ 2843.403/3 


- 2844.1667/* + 185.687500/^ 


11 


/ 


- 605.5500/2 


+ 11399.82/3 


- 23374.694/* + 5449.81250/^ 


12 


/ 


- 1288.383/2 


+ 42793.33/3 


- 161291.55/* + 88986.8125/5 - 3891.03125/'^ 



Table A2. The 12 first cumulants for e = 2. 



k 


(-l)'=fc!Xfe 






1 


/ 






2 


0.500000/ 






3 


0.307692/ 






4 


0.200000/ - 


0.050000/2 




5 


0.132231/- 


0.180995/2 




6 


0.087912/- 


0.424981/2 + 0.046429/3 




7 


0.058554/ - 


0.820875/2 + 0.387161/3 




8 


0.039024/ - 


1.420263/2 + 1.889714/3 - 


0.095141/* 


9 


0.026014/ - 


2.293572/2 + 7.080389/3 _ 


1.435478/* 


10 


0.017342/ - 


3.537360/2 ^ 22.60912/3 - 


12.021342/* + 0.349690/5 


11 


0.011561/" 


5.283592/2 + 64.91179/3- 


74.270396/* + 8.293927/^ 


12 


0.007707/ - 


7.711882/2 + 172.9145/3 - 


379.11871/* + 106.0533/5 - 2.01157/^ 
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Table A3. The 12 first cumulants for g = 3. 



k i-l)" k\xk 



1 


/ 






2 


0.333333/ 






3 


0.142857/ 






4 


0.066667/ - 


0.011111/2 




5 


0.032258/ - 


0.029762/2 




6 


0.015873/ - 


0.052596/2 + 0.003872/3 




7 


0.007874/ - 


0.077199/2 + 0.024758/3 




8 


0.003922/ - 


0.102003/2 + 0.093423/3 - 


0.002627/^ 


9 


0.001957/ - 


0.126065/2 + 0.271598/3 - 


0.032822/4 


10 


0.000978/ - 


0.148870/2 + 0.673638/3 _ 


0.223264/4 + 0.002953/^ 


11 


0.000489/ - 


0.170176/2 + 1.501653/3 - 


1.107501/4 + 0.064580/5 


12 


0.000244/ - 


0.189903/2 + 3.102167/3 - 


4.502556/4 + 0.714094/5 - 0.004752/^ 



Table A4. The 12 first cumulants for g = A. 



k 


(-l)'=fc!Xfe 






1 


/ 






2 


0.250000/ 






3 


0.081633/ 






4 


0.029412/- 


0.003676/2 




5 


0.011103/- 


0.007701/2 




6 


0.004296/ - 


0.010755/2 + 0.000597/3 




7 


0.001686/ - 


0.012572/2 + 0.003044/3 




8 


0.000667/ - 


0.013302/2 + 0.009236/3 - 


0.000136/4 


9 


0.000265/ - 


0.013214/2 + 0.021700/3 - 


0.001567/4 


10 


0.000105/ - 


0.012575/2 + 0.043644/3 - 


0.009292/4 + 0.000004/5 


11 


0.000042/ - 


0.011601/2 + 0.079054/3- 


0.039239/-* + 0.000752/5 


12 


0.000017/- 


0.010458/2 + 0.132857/3 - 


0.134113/4 + 0.010015/5 + 0.000100/*^ 
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